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Abstract 

Stimulated by recent development of the calculation methods of the running 
quark masses m q (fM) and renewal of the input data, for the purpose of making 
a standard table of rn q (/i) for convenience of particle physicists, the values 
of m q (fj,) at various energy scales //(// = 1 GeV, fj, = m c , fi = nib, M = m t 
and so on), especially at \i = mz, are systematically evaluated by using the 
mass renormalization equations and by taking into consideration a matching 
condition at the quark threshold. 
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I. INTRODUCTION 



It is very important to know reliable values of quark masses m q not only for hadron 
physicists who intend to evaluate observable quantities on the basis of an effective theory, 
but also for quark-lepton physicists who intend to build a model for quark and lepton 
unification. For such a purpose, for example, a review article [0 of 1982 by Gasser and 
Leutwyler has offered useful information on the running quark masses m g (/x) to us. However, 
during the fifteen years after the Gasser and Leutwyler's review article, there have been some 
developments in the input data and calculation methods: the QCD parameter A^r has been 
revised @; top-quark mass m t has been observed 0-0]; the three- loop diagrams have been 
evaluated for the pole mass M"^ ole || and for the running quark mass m q {p) [0]; a new 
treatment of the matching condition at the quark threshold has been proposed ||. On the 
other hand, so far, there are few articles which review masses of all quarks systematically, 
although there have been some re-estimates H~CT for specific quark masses. For recent one 
of such few works in systematical study of all quark masses, for example, see Ref. |T{| by 
Rodrigo. We will give further systematical studies on the basis of recent data and obtain a 
renewed table of the running quark mass values. 

The purpose of the present paper is to offer a useful table of the running quark masses 
m q {p) to hadron physicists and quark-lepton physicists. In Sec. [TV], by using the mass 
renormalization equation (4.1), we will evaluate the value of m q (fi) at various energy scales 
H, e.g., [A — 1 GeV, /x = m q (q = c,b,t), \i = M^ ole , \i = mz, A* = -^W, and so on, where 
M^ ole is a "pole" mass of the quark q, and A\y is the symmetry breaking energy scale of the 
electroweak gauge symmetry SU(2)i x U(l)y: 

A w = (0°) = (V2G F )'^/V2 = 174.1 GeV . (1.1) 

In the next section, we review the light quark masses m u (fi), m^(/i) and m s (n) at fi = 1 
GeV. In Sec. |T|, we review pole mass values of the heavy quark masses MP° le , Ml ole and 
M? ole . In Sec. [TV], running quark masses m q (fi) are evaluated for various energy scales /i 
below /i = A w = 174.1 GeV. In Sec. 0, we comment on the reliability of the perturbative 
calculations of the running quark masses m q (ii) (/i < Aw)- In Sec. [VI], we summarize our nu- 
merical results of the running quark mass values m g (/x), the charged lepton masses m^(/i), the 
Cabibbo-Kobayashi-Maskawa (CKM) matrix Vckm(I-i), and the SU(3) c xSU(2) L xU(l) y 
gauge coupling constants (?j(/i) {i = 1,2,3) at fi = m z . In Sec. [VI I| , for reference, the evo- 



lution of the Yukawa coupling constants is estimated at energy scales higher than /i = Ajy 
for the cases of (A) the standard model with one Higgs boson and (B) the minimal SUSY 
model. Finally, Sec. |V111 is devoted to summary and discussions. 



II. LIGHT QUARK MASSES AT // = 1 GeV 

Gasser and Leutwyler |l| have concluded in their review article of 1982 that the light 
quark masses m u (/i), m^(/i) and m s (/i) at \i = 1 GeV are 
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m u (l GeV) = 5.1 ± 1.5 MeV , 

m d (l GeV) = 8.9 ± 2.6 MeV , (2.1) 
m s (l GeV) = 175 ±55 MeV , 

from QCD sum rules. 

On 1987, Dominguez and Rafael || have re-estimated those values from QCD finite 
energy sum rules. They have obtained the same ratios of the light quark masses with those 
estimated by Gasser and Leutwyler, but they have used a new value of (m u + m d ) at /i = 1 
GeV 

(m u + m d )„ =l GeV = (15.5 ± 2.0) MeV , (2.2) 

instead of Gasser-Leutwyler's value {m u +m d ) M =i GeV = (14±3) MeV. Therefore, Dominguez 
and Rafael have concluded as 

m u (l GeV) = 5.6 ± 1.1 MeV , 

m d (l GeV) = 9.9 ± 1.1 MeV , (2.3) 
m s (l GeV) = 199 ±33 MeV . 

Recently, by simulating r-like inclusive processes for the old Das-Mathur-Okubo sum rule 
relating the e + e~ into 1 = and 1 = 1 hadron total cross-section data, Narison (1995) fK| 
has obtained the following values: 

m u (l GeV) = 4± 1 MeV , 

m d (l GeV) = 10 ± 1 MeV , (2.4) 
to s (1 GeV) = 197 ±29 MeV , 

which are roughly in agreement with (2.3). 

On the other hand, by combining various pieces of the information on the quark mass 
ratios, Leutwyler (1996) [|TTJ] has recently re-estimated the ratios 

m u /m d = 0.553 ± 0.043 , 

m s /m d = 18.9 ± 0.8 , (2.5) 

and has obtained 

m u (l GeV) = 5.1 ±0.9 MeV , 

m d (l GeV) = 9.3 ± 1.4 MeV , (2.6) 
m s (l GeV) = 175 ±25 MeV , 

The values (2.6) are in agreement with (2.1), (2.3) and (2.4). 

There is not so large discrepancy among these estimates as far as m u and m d are con- 
cerned, except for estimates by Donoghue, Holstein and Wyler (1992) JL2J, who have obtained 

m d /m u = 3A9, m s /m d = 20.7 , (2.7) 
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from the constraints of chiral symmetry treated to next-to-leading order. Eletsky and Ioffe 
(1993) and Adami, Drukarev and Ioffe (1993) |14J] have obtained 



(m d - m u ) M=0 .5 cev = 3 ± 1 MeV , (2.8) 

from the QCD sum rules on the isospin-violating effects for D and D* and for N, E and H, 
respectively. The value (2.8) is consistent with (2.3) and (2.6). The value 

(m u + m d ) M=1 GeV = (12 ± 2.5) MeV , (2.9) 

obtained from QCD finite energy sum rules and Laplace sum rules by Bijnens, Prades and 
Rafael (1995) pf is consistent with (2.2). 



On the contrary, for the strange quark mass m s , two different values, m s ~ 175 MeV, 
[(2.1) and (2.6)], and m s ~ 200 MeV , [(2.3) and (2.4)], have been reported. Recently, 
Chetyrkin et al. (1997) [|I1| have estimated 

m s (l GeV) = 205.5 ± 19.1 MeV , (2.10) 

by an order-a^ determination from the QCD sum rules. The value (2.10) is consistent with 
(2.3). (Of course, if we take their errors into consideration, these values are consistent.) 

Hereafter, as the light quark masses at \i = 1 GeV , we will use the following values 
which are weighted averages of the values (2.3), (2.4), (2.6) and (2.10). 

m u {l GeV) = 4.88 ± 0.57 MeV , 

m d (l GeV) = 9.81 ± 0.65 MeV , (2.11) 
m s (l GeV) = 195.4 ± 12.5 MeV . 



III. HEAVY QUARK MASSES 

A. Charm and bottom quark masses 

Gasser and Leutwyler (1982) |TJ have estimated charm and bottom quark masses m c and 
rub from the QCD sum rules as 

m c (m c ) = 1.27 ±0.05 GeV , 

m h (m h ) = 4.25 ± 0.10 GeV . (3.1) 

Titard and Yndurain (1994) [T7|] have re-estimated charm and bottom quark masses by 
using the three-level QCD and the full one-loop potential. They have concluded that 

MP° le = 1.570 ± 0.019 =F 0.007 GeV , 

Mf e = 4.9061°;^ =F 0.004±8SJ GeV , (3.2) 
m c (m c ) = 1.306±g-g^ ± 0.006 GeV , 

m b (m b ) = 4.397^™n^ GeV , (3.3) 
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where the first- and second-errors come from the use of the QCD parameter A^L = 0.20ig;g| 
GeV and the gluon condensate value (a s G 2 ) = 0.042 ± 0.020 GeV 4 , respectively, and the 
third error denotes a systematic error. 

On the other hand, from the QCD spectral sum rules to two-loops for ip and T, Narison 
(1995) has estimated the running quark masses 

m c (Mf T2 ) = 1.23i{j;ig ± 0.03 GeV , 

m 6 (M 6 p T2 ) = 4.23+8:81 ± 0.02 GeV , (3.4) 
corresponding to the short- distance perturbative pole masses to two-loops 

M C PT2 = 1.42 ± 0.03 GeV , 

M PT2 = 4.62 ± 0.02 GeV , (3.5) 
and the three-loop values of the short-distance pole masses 

M C PT3 = 1.64i8$ ± 0.03 GeV , 

M PT3 = 4.87 ± 0.05 ± 0.02 GeV . (3.6) 

The values (3.6) are in agreement with the values (3.2) estimated by Titard and Yndurain 
while the values (3.5) are not so. Narison asserts that one should not use M^ T3 because the 
hadronic correlators are only known to two-loop accuracy. 

Although we must keep the Narison's statement in mind, since we use the three-loop 
formula (4.5) for the running quark masses m q (fi) for all quarks q = u,d,---,t, hereafter, 
we adopt the following weighted averages of (3.2) and (3.6), 

MP° le = 1.59 ± 0.02 GeV , 

M poie = 4 g9 ± Q5 GeV ? ? ) 

as the pole mass values. 
B. Top quark mass 

The explicit value of the top quark mass was first reported by the CDF collaboration 
(1994) [[| from the data of pp collisions at \/s = 1.8 TeV: 

m t = 174 ± 10t\l GeV • (3-8) 
They (1995) [|J have also reported an updated value 

m t = 176 ± 8 ± 10 GeV . (3.9) 
On the other hand, the DO collaboration H has reported the value 

m t = 199±|? ± 22 GeV . (3.10) 
The particle data group (PDG96) |2l] has quoted the value 

m t = 180 ± 12 GeV , (3.11) 
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as the top quark mass from direct observations of top quark events. 
Hereafter, we use the value (3.11) as the pole mass of the top quark. 

C. Mass values m ? (/i) at // = MP° le 

The relation between the pole mass MP° le and the running quark mass m q (MP ole ) at 
fj, = MP ole has been calculated by Gray et al ||: 

m q (M p q ole ) = M p q ole 

where K c = 14.5, Kf, = 12.9 and K t = 11.0. The definition of K q and their estimates are 
given in Appendix A. The values of a s (fi) at various values of \x and errors are given in Table 
VII in Appendix B. By using (3.12), from (3.7) and (3.11), we obtain 

m c {MP° le ) = 1.213 ± 0.018;°;°;*° GeV , 

m b (Mf e ) = 4.248 ± 0.046:™ GeV , (3.13) 
m t (Mf ole ) = 170.1 ± 11.4 =f 0.3 GeV , 

where the first and second errors come from ±AMJ oie and ±AA|^, respectively. 



1 + 



4a s (MP° le 
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+ 0(a 3 s ) 



(3.12) 



IV. BEHAVIORS OF m 9 (/x) AT THE QUARK THRESHOLDS 



The scale dependence of a running quark mass is governed by the equation 



/^™ 9 (^) = ~l{ois)m q {n) 



(4.1) 



where 
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Then, m q (/j) is given by 



m q (n) = R(a s (fi))m, 



(4.4) 
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(4.5) 



where m q is the renormalization group invariant mass which is independent of ln(/i 2 /A 2 ), 
a s (/x) is given by (B4) in Appendix B and (i = 0, 1, 2) are also defined by (B3). By using 
(4.5) and obtained in Appendix B, we can evaluate R( n \fi) for a < fj, n +i, where \x n is 



the nth quark flavor threshold and we take \i n 

Quite recently, the four-loop beta function and quark mass anomalous dimension have 
been obtained by several authors [22]. In this paper, we evaluate the running quark masses 
by using the three-loop results (4.1)-(4.5). The effects of the four-loop results to the three- 
loop results will be discussed in the next section. 

We can evaluate the values of m q {m q ) (q = c, b, t) by using the values of M"^ ole given in 



Sec. Ill and the relation 



m q M = [r^^)/R^{M^)] m qn (M^) (/. < /i n+1 ) 



(4.6) 



Similarly, we evaluate the light quark masses m q {m q ) (q = u, d, s) by using the relation 



m q (fi) = [R {3) (fi)/R {3) (l GeV)] m q (l GeV) (ji < fi 4 ) 



(4.7) 



and the values m q (l GeV) given in (2.11). The results are summarized in Table |. The values 
of m u (m u ), m d {m d ) and m s {m s ) should not been taken rigidly, because the perturbative 
calculation is not reliable for such a region in which a s {fx) takes a large value (see the next 
section) . 



Exactly speaking, the estimates of A^ in Table |V11| in Appendix [B] are dependent on 
the choices of the quark threshold \i n = m qn (m qn ). The values in Table [VIl] and Table | 
have been obtained by iterating the evaluation of A^ and m q (m q ). 

Running quark mass values m qn (n) at \i > fi n+ i cannot be evaluated by the formula (4.4) 
straightforwardly, because of the quark threshold effects. As seen in Fig. [TJ, the behavior of 
i?(/i) is discontinuous at \i = \i n = m qn {m qn ). 

The behavior of the nth. quark mass m^} (n < N) at n N < n < fiN+i are given by the 
matching condition [BJ 



m 



(N) 
qn 



00 



m 



(JV-l) 



qn 



(m) 



1 / 2 5 89 \ 
12 [ Xn + 3 XN + 36 J 



a 



(N) 



(/*)' 



(4.8) 



where 



x N 



In (mjy^)) 2 /// 5 



(4.9) 



For example, the behavior of m c (/i) at a < /a 5 , m^(fi), can be evaluated by using (4.6), 
while those at /is < \x < ji§ and /Iq < ji must be evaluated by using (4.8) with rrv^\ji) and 

[mf\jjifj //i 2 , respectively. In 

u, d, s), where 



;r 5 



In 



KV)) Itf 



and with m^(fi) and Xq = In 



Fig. 2, we illustrate the /^-dependency of the light quark masses m q {ji) (q 
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we have taken the matching condition (4.8) into account. We can see that the discontinuity 
which was seen in Fig. [I] disappears in Fig. [| 

We also illustrate the behavior of the heavy quark masses m q (fi) (q = c,b,t) in Fig. 
H Exactly speaking, the word "the running mass value mg(/i)" of a heavy quark Q at a 
lower energy scale \i than [i = rriQ(mQ) loses the meaning. For example, the effective quark 
flavor number n q is three at fi — 1 GeV, so that the value of m t (fi) at \i — 1 GeV has 
not the meaning. However, for reference, in Fig. 3, we have calculated the value of mg(/x) 
(Q — Qn) a t A*n < A 4 < fJ'n+i ( n < -^0 by using the relation mg(yu) = ttiqR^^) [not 
m Q (//)=m Qj RW(//)]. 

The numerical results are summarized in Table [TT]. As stressed by Vermaseren et al. [22], 
the invariant mass m q is good reference mass for the accurate evolution of the MS quark 
masses to the necessary scale \x in phenomenological applications. The values of m q are also 
listed in Table II. 



V. RELIABILITY OF THE PERTURB ATIVE CALCULATION BELOW pt ~ 1 GeV 



As we noted already, the values of the light quark masses m q (m q ) (q = u,d,s) should 
not be taken rigidly, because the perturbative calculation below fi ~ 1 GeV seems to be not 
reliable. 

In order to see the reliability of the calculation of a s (fi), in Fig. f|, we illustrate the values 
of the second and third terms in { } of (B4) in Appendix [FJ separately. The values of the 
second and third terms exceed one at fi ~ 0.42 GeV and \i ~ 0.47 GeV, respectively. Also, 
in Fig. [|, we illustrate the values of the second and third terms in { } of (4.5) separately. 
The values of the second and third terms exceed one at \i ~ 0.58 GeV and /i ~ 0.53 GeV, 
respectively. These means that the perturbative calculation is not reliable below /x ~ 0.6 
GeV. Therefore, the values with asterisk in Tables I, II and VI should not be taken strictly. 

These situations are not improved even if we take the four-loop correction into consider- 
ation. For example, for n q = 3, d(a s / it) /d In fi is given by |22| 



d(a s /ir) 9 fa 



cHn/i 2 V 7r 



1 + 1.79— + 4.47 f — ) +21.0(^) + 



(5.1) 



Since the value of a s ji\ is a s /iT ~ 0.16 at /i ~ 1 GeV, the numerical values of the right-hand 
side of (5.1) becomes 

diajii) 9 /a,\ 2 r n . 

1 1 [1 + 0.28 + 0.11 + 0.085 + •••] , (5.2) 



d In /i 2 V tt 

so that the fourth term is not negligible compared with the third term. This suggests that 
the fifth term which is of the order of (a s /7r) 6 will also not be negligible below /i ~ 1 GeV. 

However, we consider that the evolution of m ? (/i) above ji ~ 1 GeV (from /i ~ 1 GeV to 
fi ~ m z ) is reliable in spite of the large error of a s (/i) at \i ~ 1 GeV. 
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VI. OBSERVABLE QUANTITIES m q (p), V C km{p) AND ai (ji) AT y, = m z 



For quark mass matrix phenomenology, values of m q (p) at fi — m z are useful, because 
the observed CKM matrix parameters \Vij\ are given at /i = m^. We summarize quark and 
charged lepton masses at \i = mz 



m u = 2.33±°;^ MeV 
m d = 469ig;|g MeV 



m e = 0.48684727 
±0.00000014 MeV , 



m c = 677±|f MeV , m t ■■ 
m s = 93.4i^;g MeV , m b 
m M = 102.75138 m T 
± 0.00033 MeV , 



= 181 ± 13 GeV , 
: 3.00 ±0.11 GeV 
= 1.74669 

+0.00030 p c v 
-0.00027 ' 



where the running charged lepton masses mi(fi) have been evaluated from the relation for 
the physical (pole) masses Mi [p3| 



a(p) 



71 



4 m 



The value of m b {mz) in (6.1) is in good agreement with the value p4 
m b (m z ) = 2.67 ± 0.25 ± 0.27 ± 0.34 GeV , 
which has recently been extracted from CERN LEP data. 



(6.2) 



(6.3) 



On the other hand, the standard expression p5[ of the CKM matrix V is given by 



V 



C12C13 



S12C13 



-S12C23 - ci 2 S23Si 3 e J<513 C12C23 - si2S 23 si 3 e i ' 513 



si 3 e 
S23C13 



\ S12S23 - C 12 C 2 3S 13 e lSls -C12S23 - Si 2 C23Si 3 e J<513 C23C13 



(6.4) 



The observed values \V US \, \V ub \ 0] an d \V cb \ |p6| , p7 



arc 



\V US \ 

\Vcb\ 

\v ub /v cb \ 



0.2205 ±0.0018 , 
0.0373 ±0.0018 , 
0.08 ±0.02 , 



(6.5) 



where the value of \Vcb\ has been obtained by combining the OPAL97 value |26j \V cb \ = 
0.0360 ± 0.0021 ± 0.0024 ± 0.0012 and the ALEPH97 value J27J | V^ fe | = 0.0344 ± 0.0016 ± 
0.0023 ± 0.0014 with the PDG96 value \V cb \ = 0.041 ± 0.003. Because of the hierarchical 
structure (V^l 2 |K&| 2 ^ |Kib| 2 , the following expression of V will also be useful: 



V 



2-A 2 
-A 



A 



-15 



1-±A 2 p 



(6.6) 



\p — oe lS —p 1 



where A = IK, 



1 = \V cb \ and a = \V ub \. Hereafter, we will use the observed values (6.5) 
as the values of |Vy(/z)| at fi = mz- Then, from the expression (6.4) (not the approximate 
expression (6.6)), we obtain the numerical expression of V{p) at [A — m z , 
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V(m z ) 



/ 0.9754 0.2205 0.0030e" i<5 

-0.2203 - 0.0001e i<5 0.9747 0.0373 j . (6.7) 

V 0.0082 - 0.0029e i<5 -0.0364 - 0.0007e i<5 0.9993 



Since we have already known the numerical values of D u = diag(m u , m c , m t ), D d = 
diag(m d ,m s ,mb) and Vij (except for the parameter 8) at p, — m z , by using the relations 

U U L M U U U J = D U1 U d L M d U% = D d , V = U2U?, (6.8) 

we can determine the numerical structures of the squared mass matrices H u and H d which 
are defined by 

H u = M U M\ , H d = M d M d . (6.9) 

Especially, at a special quark-family basis on which the up-quark mass matrix takes a diag- 
onal form D u , we can readily obtain the matrix form H u and H d : 

m 2 Jm 2 t 
H u = D 2 u = m 2 t | ml/ ml | , (6.10) 
1 



/ a 2 {l + x 2 ) pa{y + e- iS ) ae- iS 
H d = VDjV^ ~ m 2 pa(y + e iS ) p 2 (l + y 2 /x 2 ) p \ , (6.11) 
V ae iS p 1 



x = , y = — [ — . (6.12) 



where 



a mb pa \mj,> 

Numerically by using (6.7), but without using the approximate expression (6.11), we obtain 

1.66 x 10~ 10 
H u (m z ) = m 2 t (m z ) j 1.40 x 10" 5 | , (6.13) 

1 



H d {m z ) = m 2 b (m z ) 



5.84 x 10 

A5\ 



-5 



(2.08 + l.lle 16 ) x 10 
2.98 x 10- 3 e i<5 



-4 



(2.08 + l.lle- iS ) x 10' 4 2.98 x lO" 3 ^ 
2.31 x 10" 3 
3.72 x 10~ 2 



e ~\ 
3.72 x 10~ 2 
0.9986 J 



(6.14) 



where m 2 (m z ) = 3.24 x 10 4 GeV 2 and m 2 b (m z ) = 9.00 GeV 2 . 

In the standard model [not SU(2) L x SU(2) R x U(1) Y , but SU(2) L x U(1) Y ], by a suit- 
able transformation of the right-handed fields, we can always make quark mass matrices 
(M u , M d ) Hermitian. Furthermore, in the quark-family basis where M u = D u , the quark 
mass matrices are given by 
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M u = D u = m t 



I m u /m t 
m c /m t 
V 1 



(6.15) 



M d = VD d V^ 



m s (m d 



m b \m s 
m 



m 



+ \ A ) A— ae 



X- 



V 



ae 



m b 
id 



P 



P 
1 



/ 



It is well known that if we assume (M d )u = 0, we obtain the relation p8 



A = \V US \ ~ \J-m d /m s 
Then, we obtain a simpler expression of M d 

( 



M d ~ m b 



m d m s g 
s— ae 10 



mi 



m d m s m s 

— P 



mi 



V 



ae 



iS 



m b 
P 



Numerically, by using Eq. (6.7), we obtain 



1.29 x 1(T 5 
M u (m z ) = m t (m z ) | 3.75 x 10~ 3 

1 



(6.16) 



(6.17) 



(6.18) 



(6.19) 



M d (m z ) = m b (m z )x 

( 3.01 x 10~ 3 (6.36 + O.lle-^) x HT 3 (-0.24 + 2.97e~ i5 ) x 10~ 3 \ 

(6.36 + 0.11e i5 ) x 10~ 3 0.0310 0.0362' ,(6.20) 

V (-0.24 + 2.97e i<5 ) x 10~ 3 0.0362 0.9986 / 

where m t (m z ) = 180 GeV and m b (m z ) = 3.00 GeV. For the case of m s < 0, instead of 
(6.20), we obtain 

M d (m z ) = m b (m z )x 

I -1.8 x 10~ 5 (-7.03 + 0.11e~ il5 ) x 10~ 3 (0.26 + 2.98e- i<5 ) x 10~ 3 \ 

(-7.03 + 0.1 le^) x 10~ 3 -0.0281 0.0384 .(6.21) 

V (0.26 + 2.98e i5 ) x 10~ 3 0.0384 0.9986 / 
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We can obtain quark mass matrix forms on arbitrary quark-family basis by the unitary 
transformation H' u = UH U W and H' d = UH d W for (6.10) and (6. 11), respectively [and also 
M' u = UM U W and M' d = UM d W for (6.15) and (6.16), respectively]. Explicit mass matrix 
forms on other special quark-family basis are, for example, given in Refs. p9| , |30|| . 

By starting from the numerical expressions of the mass matrices H u and Hd at fi = 
rriz, (6.10) and (6.11), we can also obtain the mass matrix form M q (q = u,d) (in other 
words, the Yukawa coupling constants) at arbitrary energy scale fi which is larger than 
the electroweak scale A^. In the next section, we discuss the evolution of the Yukawa 
coupling constants. Then, we will use the following values of the SU(3) c xSU(2) L xU(l)y 
gauge coupling constants at /i = m z : 

ai(m z ) = 0.016829 ± 0.000017 , 

a 2 (m z ) = 0.0334931°;™ , (6.22) 
a 3 (m z ) = 0.118 ±0.003 . 



which are derived from [31 



a(m z ) = (128.89 ± 0.09)" 1 , 

sin 2 6 W = 0.23165 ± 0.000024 , (6.23) 

and Ag^ = 209±i§ MeV §. Here, the coupling constants of U(l)y, SU(2) L , and SU(3) C 
gauge bosons, g%, gi and g$, are defined as they satisfy the relation 

? = §2 + 7, ' (6 - 24) 

and the relation in the SU(5)-GUT [32] limit 

9i = 92 = #3 • (6.25) 



VII. EVOLUTION OF YUKAWA COUPLING CONSTANTS 



So far, we have evaluated values of the running quark masses m q (fi) at energy scales 
which are below the electroweak symmetry breaking energy scale Ky/ by using the formula 
(4.1). However, for the quark masses at an extremely high energy scale far from A^, we 
must use "evolution" equations of Yukawa coupling constants yf, (a = u,d : i,j = 1,2,3). 
The numerical results of the Yukawa coupling constants have already been given in many 
literatures. Since our interest in the present paper is in the updated values of the quark 
masses m g (/x) (i.e., the Yukawa coupling constants y q ), we give only a short review of the 
evolution of the Yukawa coupling constants, and do not give a systematical study of the 
numerical results. 

We define the Yukawa coupling constants yfj as follows: 

3 3 

H mass = Y,Y,y?$Lai^Raj(Pa + h.C. , (7.1) 

a=u,d i=l j=l 
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where (f> a are the vacuum expectation values of the neutral components of the Higgs bosons 
<j) a which couple with fermions ip a , and they mean <fi° u and <p d for the minimal SUSY model 
(Model B), while they mean a single Higgs boson <fi° u = <fi d = 0° for the standard model with 
one Higgs boson (Model A). The quark mass matrices M u and M d at /i = A w are given by 



1 

V2 



Y a (/l)v a 



(7.2) 



where Y a denotes a matrix (Y a )ij = y^, and v a are the vacuum expectation values of 



v a = y/2(<j)l), and v u = v d = V2A W for Model A and ^v 2 u + vj = V2A W for Model B. 

The evolution of the coupling constants Y a (/j,) from Y a (A w ) is given by the following 



equations [|33f| : 



dYa 

dt 



16ir 



ift 



(i) 



167T 



0, 



(2) 



Y„ 



u, d, e) 



(7.3) 



t = \n(fx/A w ) 



b 



(7.4) 
(7.5) 



ft 



(2) 



b.c 



(7.6) 



=yy t 



(7.7) 



where, for convenience, we have changed the definition of the Hermitian matrix H a from (6.8) 
in the previous section to (7.7). The coefficients c a 1 ^ and a b a in the one-loop contributions 
(3^ are given in Table |TTT| according to Models A and B, where 



,(i) 



(7.8) 



The coefficients c a 2 \ b b a and b b a c in the two- loop contributions (3^ are given in Appendix C, 
because they have too long expressions. The evolution of the gauge coupling constants g^fi) 
is given in Appendix |D| 

By using the information of Vij(fi) at fi — mz in the previous section, we can obtain not 
the knowledge of M q (mz), but that of H q (mz), i.e., H u = D\ and H d = VD d V^ ', where D q 
(q = u, d) are the diagonalized matrices of Y q . Then, the expression for H a (/x) 



dt 



-H„ 



16n' 



(167T 2 ) 2 



ft 



(2) 



H a + H a 



16n 



0, 



(i)t 



2^a 



+ 



(16tt 



2^2 



ft 



(2)t 



(7.9) 



is useful rather than (7.3) which is the expression for Y a . Hereafter, for simplicity, we 
calculate the evolution not from \i = A w , but from \i = m z because most of the input 
values at /i = mz have already given in Sec. [VI]. Since the numerical results are insensitive 
to the value of the phase parameter 5 13 (tt/3 < 5^ < 2tt/3) in the CKM matrix V, (6.4), 
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we will use the value 5 = 5 is = tt/2 below. For Model A (Standard model with one Higgs 
boson), we must assume the value of the Higgs boson mass m H . We will take a typical value 
m H = V2A W = 246.2 GeV (see later discussion). For Model B (Minimal SUSY model), 
we must assume the value of tan/5 = v u /vd- We will take a typical value tan/3 = 10. The 
numerical results of y q are given below. Here, the values y^ are obtained by diagonalizing 



the matrix H a , and it does not mean J (H a )u- 

(A) Standard model with one Higgs boson 

As seen in Appendix |A], in the calculation of the two-loop contributions, the evolution of 
the Yukawa coupling constants y q depends on the coupling constant Xh of the Higgs boson 
4>, which is related to the Higgs boson mass mn as 

\ H = m 2 H /v 2 . (7.10) 

We find [53 1 that the input value of mj/(m^) which is less than 2.2 x 10 2 GeV leads to a 



negative Xh at a unification scale /x = M x , while that which is larger than 2.6 x 10 2 GeV 
leads to the burst of Xh at the unification scale. Therefore, if we put an ansatz that Nature 
accepts only the parameter regions in which the perturbative calculations are valid, we can 
conclude that the Higgs boson mass vhh in the standard model must be in 

220 GeV < m H (m z ) < 260 GeV . (7.11) 

In Table [TV], we list the numerical results of m q (/i) = y q {\i)vj \/2 at the typical energy scales 
H = m z , n = 10 9 GeV and \i = M x . For the comparison with the SUSY model (Model 
B) later, the values m q (fi) at \i = M x are listed, where M x is a unification scale of SUSY, 
M x = 2 x 10 16 GeV. Here, we have tentatively taken a value m H = y/2A w = 246.2 GeV 
(i.e., Xh = 1) as the input value of mj(mz). 

We also obtain the numerical expression of the CKM matrix V(fi) at \i — M x 

( 0.9754 0.2206 -0.0035A 

V(M X ) = -0.2203 0.9745 0.0433 , (7.12) 

\0.0101e- 19Oi -0.0422e +10 °' 0.9991 / 

correspondingly to (6.7) at fi — mz, where we have taken 5 = 90° tentatively. We also 
obtain the numerical result of (M u ,M d ) at \i = M x correspondingly to (6.19), (6.20) and 
(6.21): 



M U (M X ) = m t {M x ) 3.23 x 10~ 3 , (7.13) 




M d {M x ) = m h {M_ 



x, 



I 0.0035 

0.0074e +1 - 2Oi 
V 0.0035e +95 - 3Oi 



0.0074e- L2Oi 

0.0363 
0.0418e-° 03Oi 



0.0035e- 95 ' 3Oi 
0.0418e +a03Oi 
0.9982 



(7.14) 



and 
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-1.9 x 1(T 5 -0.0082e +11 °* 0.0035e" 84 - 1Oi 
M d (M x ) = m b (M x ) | -0.0082e- L1Oi -0.0324 0.0447e- - 040 * | . (7.15) 

0.0035e +84 - 1Oi 0.0447e +004Oi 0.9980 

where m t (M x ) = 84.2 GeV and m b (M x ) = 1.071 GeV. 

(B) Minimal SUSY model 

The scale of the SUSY symmetry breaking msusY is usually taken as msusY — m t or 
msusY — 1 TeV. For simplicity, we take msusY = °mz in the present numerical study, 
because the numerical results of y g (fi) are not sensitive to the value of msusY- 

The values of m q (fi) = y g (fi)v/^/2 (q = u, d) are sensitive to the value of tan/3 = v u jv d . 
A large value of tan/?, tan/3 ~ 60, leads to the burst of m(,(/i) at the unification scale 
fi = M x — 2 x 10 16 GeV. On the other hand, a small value of tan/3, tan/3 ~ 1.5, leads to 
the burst of m t {fi) at the unification scale. The values of m q (fi) are insensitive to the value 
of tan /3 in the region from tan (3 ~ 5 to tan (3 ~ 30 [ 35| . In Table |V|, we list the numerical 



results of m q (fi) at the typical energy scales, fi = mz, fi = 10 9 GeV and fi = Mx- Here, we 
have tentatively taken a value tan (3 = 10 as the input value of tan f3. 

In Fig. |6|, for reference, we illustrate the behavior of mt(n), mb(fi) and m T (fi). The 
value of m t (Mx) is highly dependent on the input value of m t (mz)- Therefore, the value 
of m t (M x ) in Table [V] should not be taken strictly. Also, the energy scale fix at which 
mbi^x) = ' m r(/ i x) is highly dependent on the input value of m h {m z ). Therefore, the value 
fix should also not be taken strictly. 

As seen in Fig. ^ it is very interesting that the observed top quark mass value is given 
by almost the upper value which gives m q (Aw) < rn(M x )- However, since the purpose of 
the present paper is not to investigate the evolution of the Yukawa coupling constants in 
the SUSY model under some postulation [e.g., m&(//) = m T (fi) at fi = Mx], we do not go 
further more. Some of such studies will be found in Refs. p5|j36| . 

We also obtain the numerical expression of the CKM matrix V(fi) at fi = Mx 

I 0.9754 0.2205 -0.0026A 

V(M X ) = -0.2203e+°- 03Oi 0.9749 0.0318 , (7.16) 
V 0.0075e~ 19Oi -0.0311e +1 ° Oi 0.9995 / 

correspondingly to (6.7) at fi — mz-, where we have taken 5 = 90° tentatively. We also 
obtain the numerical result of (M u ,M d ) at fi = M x correspondingly to (6.19), (6.20) and 
(6.21): 



M U (M X ) = m t (M x ) 2.33 xl0- 3 , (7.17) 





/ 0.0026 0.0054e- a9Oi 0.0025e~ 93 - 9Oi \ 

M d (M x ) = m b (M x ) 0.0054e +a9Oi 0.0263 0.0310e +a03Oi , (7.18) 

V 0.0025e +93 9Oi 0.0310e- a03Oi 0.9990 / 

and 
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-1.6 x 1(T 5 -0.0060e +0 - 8 °* 0.0026e" 85 - 8Oi 

M d (M x ) = m b (M x ) | -0.0060e-°' 8Oi -0.0241 0.0326e-° 030 * | . (7.19) 

0.0026e +85 - 8 °* 0.0326e +003Oi 0.9990 

where m t (M x ) = 129.3 GeV and m b (M x ) = 0.997 GeV. 



VIII. SUMMARY 



In conclusion, we have evaluate the running quark mass values m q (n) (q = u,d, s, c, b, t) 
at various energy scales //(// = 1 GeV, /i = m q , /i = m z , and so on). The values of m q {m q ) 
given in Table |I| in Sec. |V] will be convenient for hadron physicists who want to calculate 
hadronic matrix elements on the bases of quark-parton model, heavy-quark effective theory, 
and so on. Also, the values of m g (/i), m^/i), |\4,-(/i)| and ct^/i) at \i = m z given in Sec. |Vj 
will be convenient for quark and lepton mass-matrix model-builders. In quark mass matrix 
phenomenology, the values of m q {ii) at /j, — 1 GeV have conventionally been used. However, 
we recommend the use of the values m q {m z ) rather than m q {l GeV), because we can use 
the observed values of \Vij\ as the values |^-(m^)| straightforwardly, and, exactly speaking, 
the value of m t (l GeV) does not have the meaning. 

Although, in Sec. |V11| , we have given the values of m q (fi) at n = Mx, i.e., the evolution 
of the Yukawa coupling constants y q (fi), the study was not systematical in contrast to the 
study for /i < A w . The values of y q (fi) in the standard model with one Higgs depend on the 
input value of the boson mass m^(m^). The values of y q (fi) in the minimal SUSY model 
depend on the values of the parameters msusY and tan/3. Therefore, the values m q (Mx) 
given in Table [IV] and Table [V] in Sec. |V1J| should be taken only for reference. 

We hope that the most of the present results, Table |lj in Sec. [TV] and (6.1), (6.7), (6.13) 
and (6.14) in Sec. jVJ are usefully made by particle physicists. 
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APPENDIX A: RELATION BETWEEN m q (m q ) AND M? oie 

The pole mass, M pole (p 2 = m 2 ), is a gauge-invariant, infrared- finite, renormalization- 
scheme-independent quantity. Generally, mass function M(p 2 ), which is defined by |lj 



S(p) = Z(p i )/(M(p 



(Al) 



Z(p 2 ) = 1 - %L ( a - 36 + ~ ) X H + 0(a 2 : 



(A2) 



is related to 



M(p 2 ) = m(ji) 



1 + — (a + A6) + 0(a 2 s 

7T 



(A3) 



4 m 2 m 2 — p 2 m 2 — p 2 
a = - - In — H - In — 

3 jj z p z mr 



(A4) 



m 2 — p 2 / m 2 m 2 — p 2 " 
1 + — In — 



3p 2 



(A5) 



where Xh is given by Xh = in the Landau gauge and by Xh = 1 in the Feynman gauge. 
For p 2 = m 2 , we obtain a = 4/3 and b = 0, so that we obtain the relation 



M? ole (p : 



,. (/ w . m~) = m q {m q ) [ 1 + --^ + 0(a 2 s ] 



2 2\ 



4 a, 



(A6) 



Similarly, for the spacelike value of p 2 , p 2 = — m 2 , we obtain a = 4/3 — 2 In 2 and b 
(2/3)(l — In 2), so that we obtain the gauge- dependent "Euclidean" masses 



Mf > 2 = -m 2 ) = m q (m q ) 



1 + f g-21n2)+0(«r 



(A7) 



The estimate of the pole mass has been given by Gray et al (also see 



m 



q {Mf e ) = M p q olt 



4 a s (MP° le 

+ 3 7T 



+ K a 



a 



( M pole 



IT 



+ 0(a 3 s ) 



(A8) 



A n ~ 1 

6 i=i 



(A9) 
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1 2, « 7 2 1>/„n 3673 / 1 2 71 

Kn = -7T 2 In 2 + — 7T 2 - -C(3) + — 7T 2 + ) 71 . 

9 18 6 SV ; 288 V18 144/ 



(AlO) 



A(r) 



In 2 r H — 7T 2 — fin r H — ) r 2 
6 V 2/ 



-(1 + r)(l + r 3 )L + (r) - (1 - r)(l - r 3 )L_(r) 



(All) 
(A12) 



L±(r) 



o x ± 1 



(A13) 



Here the sum in (A9) is taken over n — 1 light quarks with masses Mf°' e (Mf ole < MP° le = 



MP° le ). The numerical results are summarized in Table [VI 



In Table |VT], the values of MP° le and m q (MP° le ) for the light quarks q = u,d,s have been 
obtained by solving the relations (A8) with the help of (A7) with the inputs (2.11). These 
values for the light quarks should not be taken rigidly, because the perturbative calculation 
is unreliable for the region at which a s (n) takes a large value. Fortunately, the values of K q 
are not sensitive to the values of M^ ole for the light quarks q — u,d, s. Therefore, the values 



of K q in Table [VI| are valid not only for the heavy quarks q = c,b,t but also for the light 
quarks q — u,d, s. 



APPENDIX B: ESTIMATE OF A 



MS 



The effective QCD coupling a s = g\j Att is governed by the /3-function: 

da.. 



where 



/'■ 



at \ @° 2 



dfj, 



(3(a s 



47T 



2 a s 



64vr 3 



at + 0{a\ 



(Bl) 



(B2) 



2 19 5033 325 

0o = 11 - -n q , ft = 51 - — n„ , ft = 2857 - — — 



-n a + 



3 q ' Mi 3 9 ' r ' 9 9 ' 27 '" 9 ' 



(B3) 



and n q is the effective number of quark flavors [39]. The solution a s (fi) of (Bl) is given by 
[2] 



. 4*1 J 2ft InL 4/Jf 



2/ 8ft 2 4 



+ 



L 3 



(B4) 



where 
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L = ln(/i 2 /A 2 



(B5) 



The value of a s (fj,) is not continuous at nth quark threshold [i n (at which the nth quark flavor 
channel is opened), because the coefficients /3 , Pi and f3 2 in (B2) depend on the effective 
quark flavor number n q . Therefore, we use the expression a^(fi) (B3) with a different 
for each energy scale range fi n < \i < fi n +i- The relationship between A^ and A^- is 



fixed at /i = where m q n ^ is the value of the nth running quark mass 
and is given as follows [40]: 



in 



(») 



^qny^qn ) > 



(n-l) 



In 



A 



(r0 
MS 



(n-l) 



A. 



(n) /3(n-l) 

Po 



+2 



MS 

jfr) A (n-n 

?(«) 



In L 



■(n) 



MS 



.(n-l)' 



2$ 

A) 



V MS 
(n-l) 



(n-l) 



hi 



(n) 



(n) 



/o(«-l) 

Po 



2/3) 



(n-l)' 



^MS 
2 



A 



(n) 



(n-l) 



A> 



+ 



,A 



(n-l) 



A 



(n-l) 



(n-l) 
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-(n) 
V MS 



(B6) 



where 



rn 



(«) 



ms7 



(B7) 



Particle data group (PDG96) [2] has concluded that the world average of A^ is 



A 



<5L 

MS 



209j|jMeV 



(B8) 



Starting from A^ = 0.209 GeV, by using the relation (B6), at /15 = m^m^) = 4.339 GeV, 



(n) 

170.8 GeV, we evaluate the values of A^ 



fi4 = m c (m c ) = 1.302 GeV, and = m t (m t ) 
for n = 3, 4 and 6. The results are summarized in Table |V11| . 

We show the threshold behaviors of qW(/j) in Fig. |7|. We can see that 
A*n-i < A* < £*n connects with a™(|u) in //„ < /i < /i n+ i continuously. 



a 



(n-l) 



(//) in 



APPENDIX C: EVOLUTION OF THE YUKAWA COUPLING CONSTANTS 



The coefficients c£\ b b a and b h £ in the two-loop contributions (3^ are given as follows. 



Here, T a (a = it, d, e) are given in Table [nj in Sec. |VII| and n g is the number of generations 



(A) Standard model with one Higgs scalar 



idd 



3 
2 



11 

T 



3 
2 



!-,(/ U 



4 ' u « 



iud 



(CI) 
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K = ~\t„ - 6A H + + if A + 16 9 f , 

bt = - j2i " 6A H + ^g? + if s? + 16gf , 

63 = *7i - 2A H - - ^ + 16 9 |) , 

^ = + ^ + §AW (so 

19 2 o /35 \ 4 22 /404 80 \ 4 

+ 15^103 - ^ - %J + 9#2#3 - - J #3 , 

c f = _x 4 + -y 4 + -A 2 , - f — + —n g ) g\ - —(fig* 
d 2 2 V 200 45 V 1 20 yiy2 

31 2 2 /35 \ 4 „ 2 2 /404 80 \ 4 

+ Y^3 - ^ - %J ^2 + 90203 - - J f3 , 

v 5^ 3 x2 / 51 11 \ 4 27 2 2 /35 \ 

X 4 + -F 4 + -Ai + 1 n„ at H gigl n 

2 2 H V200 5 V yi 20 yiy2 V 4 V 

X 4 = ^Tr (3H 2 U - 2 -H u H d + 3H 2 + , 
n = + ^ + 8^) M« + Qff? + \ 9 l + 8gi) TrH d 

+ i^/l + ^2 2 ) M e , 

= m 2 H /v 2 , 
20 



and the evolutions of gi (i = 1, 2, 3) and Xh are given in Sec. [D]. 
(B) Minimal SUSY model 



1UU idd _A idd inu _ 2 Ifie A 

u u u d ^ ' u u u d ^ 7 u e ^7 



bT = bf = -2, 6f = 6f = 0, (C8) 



b u u = -3T U + 2 -g\ + 6g 2 2 , b d d = -3T d + tg\ + 6g 2 2 , &« = -3T e + 6^ 2 2 



= -T d + -</? , 6« = -T u + -^ 2 , (C9) 



= -3Tr (3tf 2 + H u H d ) + (^ 2 + 16^ 2 ) TrH u 



/403 26 ^ 4 , ^ 2 , 136 ^ 2 



+ 1 ™« 0i + SRSb-\ ofoo 

V450 15 V yi yiy2 45 yiya 



21 R \ 4^0 2 2 /304 32 \ 4 
y - 6n s J # 2 + 8# 2 £ 3 - ^— - y n g J # 3 



= -3Tr (3#J + + tf 2 ) + (--gf + 16^) M d + -gfTrH. + (- + -n fl ) 

2 2 8 2 2 f 21 r \ 4 o 2 2 / 304 32 \ 4 
+0102 + g0103 - - Qn 9 J 92+ 80203 - y— - ) & ' 

cf = -3Tr (3H 2 d + H u H d + # e 2 ) + (~g* + 16gl) TrH d + \glTzH e 



.'27 18 \ , 9 „ „ /21 

( + y ™ 9 J 0i 4 + g0?0 2 2 - ( y - K ) 2 4 • ( < i o ) 
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APPENDIX D: EVOLUTION OF THE GAUGE COUPLING CONSTANTS 



Evolution of gauge coupling constants is given by 



dt 



16tt 2 



16tt 2 ) 2 (16tt 



2\2 



c ia TiH a 



(Dl) 



where the coefficients b{, bit- and Q a are given in Table |V11J| . 

The evolution of the coupling constants A# given in (C7) is given by 



dX 



H 



-A 



(i) 



dt " 16n 2hx (16tt 2 ) 2 



;d2) 



W _ i o \ 2 ^ 9 „2 , n \ , ^ / ^ „4 , 2 _o 2 . . i 



/3r = l2y H - ( + 9gi) \ H + i l-^gf + -g(gi + g\ 



+ AX H Tr(3H u + 3H d + H e ) - 4Tr(3# 2 + 3H 2 + H, 



e I i 



(D3) 



/3f = -78A| + 18 (jfl? + 3^ ) 



313 _ \ 4 117 2 2 9 /229 50 

10n o g% g 2 g\ H 

8 9 J 32 20 12 25 V 24 9 



!T n 9 9i 



/497 \ 6 3/97 8 \ 2 4 9 /239 40 , , , , , „ 

+ VY ~ g ) 92 ~ 5 V24 3™ 9 J M ~~ 25 l"24 + M ~~ ToF I o7 + "77^9 ) 9\ 



27 /59 40 
125 V24 + ~9 



-64^ 2 Tr (# 2 + Hi) - -glTi (2H 2 U - H 2 d + 3H 2 ) - -g 4 2 Ti{3H u + 3H d + H e 



+10A 



H 



Y 9 2 + ~/ 2 + Sg 2 ) TrH u + (irf + \g\ + 8g 2 ) TrH d + \{g\ + g 2 2 )TrH e 



+ 5^ 



~g 2 + 2lg^) TrH u + (^g\ + 9<? 2 2 ) TrH d + (-y^ 2 + ll^ 2 ) M, 



-24A^Tr(3# n + 3H d + H e ) - X H Tx 3{H U - H d ) 2 + H, 



+ 20Tr(3# u 3 + 3H 3 d + H 3 e ) - 12Tr[H u H d {H u + H d )\ . 



(D4) 
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FIG. 1. Threshold behavior of R^ n \n) versus fi. 
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FIG 2. Light quark masses m ? ( / u) (q = u, d, s) versus fi. 
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FIG. 3. Heavy quark masses m q (fi) (q = c,b,t) versus [i. 




FIG. 4. Reliability of the perturbative calculation of a s (/u). The curves show the behaviors 
of the second and third terms in { } in (B4). 
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FIG. 5. Reliability of the perturbative calculation of m q (n). The curves show the behaviors of 
the second and third terms in { } in (4.5). 




!U 1Q 4 1Q 10 1Q 16 

Energy Scale (GeV) 

FIG. 6. Behavior of the Yukawa coupling constants yt(^), Vb{p) and y T (/J.) of in the minimal 
SUSY model. For convenience, the values are illustrated by the form mt(/J.) = yt(//)usin/3/v / 2, 
mb(/J,) = Ubi^v cos /3/\/2 and m T {ji) = y T (/j,)v cos (3/ y/2. 



27 



T 




m = m c (m c ) He = m t (m t ) 

10° 10 1 10 2 10 3 

\i (GeV) 



FIG. 7. Threshold behavior of af > \n) versus fi. 
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TABLES 



TABLE I. Running quark mass values m 9 (/x) at fi = m q . Input values m q (l GeV) for q = u,d,s 
and m q (MP° le ) for q = c,b,t are used. The first and second errors come from ±Am q (or ±AMP ole ) 

and LAA^L respectively. The values with asterisk should not be taken rigidly, because these 
values have been calculated in the region with a large a s (fi). 





Input : m g (l GeV) or m q (MP° le ) 


Output : m q {m q ) 


u 

d 

s 


4.88 ± 0.57 MeV 
9.81 ±0.65 MeV 
195.4 ± 12.5 MeV 


*0 4Sfi +u - UUi +u - 058 GeV 

*U. / ±OO_ 002 _o.052 ^ ev 

*0.448 ± 0.001±8;8|§ GeV 
*0.553 ± 0.0051^1 GeV 


c 

b 
t 


1.213 ± 0.018J;^4 GeV 
4.248 ± 0.046;H36 GeV 
170.1 ± 11.4 =f 6.3 GeV 


1.302 ± 0.018J;^ GeV 
4.339 ± 0.046+8;81? GeV 
170.8 ± 11.5 =F 0.2 GeV 



TABLE II. Running quark masses m q (p) and invariant masses m q (in unit of GeV). The values 
with asterisk should not be taken strictly, because the perturbative calculation is not reliable in 
the region with a large a s (fi). 



Q = 


u 


d 


s 


C 


b 


t 


M pole 


*0.501 


*0.517 


*0.687 


1.59 


4.89 


180 




+0.068 
-0.061 


+0.068 
-0.062 


+0.074 
-0.067 


±0.02 


±0.05 


±12 


m q (MP° le ) 


*0.0307 


*0.0445 


*0.283 


1.213 


4.248 


170 


+0.0022 
-0.0026 


+0.0018 
-0.0023 


+0.013 
-0.016 


+0.052 
-0.058 


+0.082 
-0.086 


±12 


m q (m q ) 


*0.436 


*0.448 


*0.553 


1.302 


4.339 


171 


+0.059 
-0.054 


+0.060 
-0.054 


+0.064 
-0.057 


+0.037 
-0.038 


+0.073 
-0.076 


±12 


m q {\GeV) 


0.00488 


0.00981 


0.1954 


1.506 


7.18 


475 


±0.00057 


±0.00065 


±0.0125 


+0.048 
-0.037 


+0.59 
-0.44 


+86 
-71 


m q (m c ) 


0.00418 


0.00840 


0.1672 


1.302 


6.12 


399 


m c = 1.302 


+0.00056 
-0.00060 


+0.00071 
-0.00077 


+0.0137 
-0.0150 


+0.037 
-0.038 


+0.32 
-0.25 


+58 
-51 


m q (m b ) 


0.00317 


0.00637 


0.1268 


0.949 


4.34 


272 


m b = 4.339 


+0.00052 
-0.00056 


+0.00073 
-0.00081 


+0.0142 
-0.0159 


+0.063 
-0.070 


+0.07 
-0.08 


+26 

-25 


m q (m w ) 


0.00235 


0.00473 


0.0942 


0.684 


3.03 


183 


mw = 80.33 


+0.00042 
-0.00045 


+0.00061 
-0.00067 


+0.0119 
-0.0131 


+0.056 
-0.061 


±0.11 


±13 


m q (m z ) 


0.00233 


0.00469 


0.0934 


0.677 


3.00 


181 


m z = 91.187 


+0.00042 
L -0.00045 


+0.00060 
-0.00066 


+0.0118 
-0.0130 


+0.056 
-0.061 


±0.11 


±13 


m q (m t ) 


0.00223 


0.00449 


0.0894 


0.646 


2.85 


171 


m t = 170.8 


+0.00040 
-0.00043 


+0.00058 
-0.00064 


+0.0114 
-0.0125 


+0.054 
-0.059 


±0.11 


±12 


m q (A w ) 


0.00223 


0.00448 


0.0893 


0.645 


2.84 


171 


A w = 174.1 


+0.00040 
-0.00043 


+0.00058 
-0.00064 


+0.0114 
-0.0125 


+0.054 
-0.059 


±0.11 


±13 


m q 


0.00496 


0.00998 


0.199 


1.59 


7.87 


546 




+0.00095 
-0.00101 


+0.00141 
-0.00153 


+0.028 
-0.030 


+0.15 
-0.16 


+0.40 
-0.41 


±49 
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TABLE III. Coefficients /% in the evolution equations of Yukawa coupling constants Y a . 



ModelA 
StandardsingleHiggs 


ModelB 
SUSY 


Gu = %gi + \gi + 8gi 
G d = \g\ + \g\ + Sgi 
Ge = \g\ + \g\ 


Gu = %gi + Sgi + fgi 
G d = y\ + Sgi + fg 2 3 

s~i Q 9 i o 9 

c e = yt + 3^1 


T u = T d = T e 
= 3Tr(H u + H d ) + TrH e 


T M = 3TrH u 
T d = T e 
= 3TrH d + TrH e 


< = < = +3/2 

< = o3 = -3/2 
a* = +3/2 


< = <4 = +3 

< = a u d = +1 
a* = +3 



TABLE IV. Evolution of the Yukawa coupling constants y a in the standard model with one 
Higgs boson (Model A). For convenience, instead of y a (^), the values of m a (fi) = y a (fi)v/V^ are 
listed, where v = y/2Aw = 246.2 GeV. The errors ±Am at fi = 10 9 GeV and fi = mx denote only 
those from ±Am at fx = mz- 





fi = m z 


H = 10 9 GeV 


li = M x 


m u (n) 


2.331^ MeV 


1.281^ MeV 


0.94l££ MeV 


m c (n) 


677+^ MeV 


3711^ MeV 


2721^ MeV 


m t (iJ,) 


181 ± 13 GeV 


1091$ GeV 


841$ GeV 


m d (fj) 


4.691^6 MeV 


2.60lft$ MeV 


1.941&28 MeV 


m s (fx) 


93.4 + {^ MeV 


51.91^ MeV 


38.71^ MeV 


m b {n) 


3.00 + 0.11 GeV 


1.511^ GeV 


1.07 ±0.04 GeV 


m e (n) 


0.48684727 MeV 
±0.00000014 


0.51541746 MeV 
±0.00000015 


0.49348567 MeV 
±0.00000014 


™m 


102.75138 MeV 
±0.00033 


108.78126 MeV 
±0.00035 


104.15246 MeV 
±0.00033 


m T (jj) 


1746.7 + 0.3 MeV 


1849.2 ±0.3 MeV 


1770.6 ±0.3 MeV 
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TABLE V. Evolution of the Yukawa coupling constants y a in the minimal SUSY model (Model 
B). For convenience, instead of y a {ji), the values of m a (fi) = y a (/j,)v sin/3/\/2 for up-quark sector 
and m a ( / u) = y a (n)v cos P/y/2 for down-quark sector are listed, where v = The errors 

±Am at n = 10 9 GeV and \i = Mx denote only those from ±Am at \i = mz- 





/j, — mz 


H = 10 9 GeV 


H = M X 


m u (fj) 


2.331^5 M eV 


lA7t u f 8 MeV 


1.041^ MeV 


m c (/i) 


677+^ MeV 


4271^ MeV 


302+^ MeV 


m t (n) 


181 ± 13 GeV 


1491^ GeV 


1291^ GeV 


m d (n) 


4.691^ MeV 


2.281^ MeV 


MeV 


m s (ii) 


93.4i^;g MeV 


45.31^ MeV 


26.5lj$ MeV 


m b (n) 


3.00 ±0.11 GeV 


1.60 ±0.06 GeV 


1.00 ±0.04 GeV 


m e (n) 


0.48684727 MeV 
±0.00000014 


0.40850306 MeV 
±0.00000012 


0.32502032 MeV 
±0.00000009 


™M 


102.75138 MeV 
±0.00033 


86.21727 MeV 
±0.00028 


68.59813 MeV 
±0.00022 


m r (/i) 


1746.7 ±0.3 MeV 


1469.51^ MeV 


1171.4 ±0.2 MeV 



TABLE VI. Pole masses M^ ole and the related quantities. The values with asterisk should not 
be taken rigidly, because these values have been calculated in the region with a large a s {jj). 





K 


A(M,/M n ) 


K 


M pole 


m q (MP° le ) 


u 


16.11 





♦ 16.11 


♦0.501MeV 


♦0.0307MeV 


d 


15.07 


♦0.838 


♦16.19 


♦0.517 MeV 


♦0.0445 MeV 


s 


14.03 


♦ 1.364 


*15.85 


♦0.687 MeV 


♦0.283 MeV 


c 


12.99 


1.114 


14.47 


1.59 GeV 


1.213 GeV 


b 


11.94 


0.746 


12.94 


4.89 GeV 


4.248 GeV 


t 


10.90 


0.0555 


10.98 


180 GeV 


170.1 GeV 



TABLE VII. The values of A-^ in unit of GeV and a s (fi n ). The underlined values are input 
values. 



n 


A^ 






3 


n QQQ+0.U4Y 
U.OOO_ .042 


1.69^ 


^3 = 0.553 GeV 


4 


0.291™ 


0-379±8:^ 


Ha = 1.302 GeV 


5 


0.209^ 


n OOO+0.013 
u.zzz_ 012 


/x 5 = 4.339 GeV 


6 


0.0882™ 


0.10 /8_ 0035 


Us = 170.8 GeV 
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TABLE VIII. Coefficients in the evolution equations of gauge coupling constants. 



Model(A) 


Model(B) 


fe i = " (to + K) 

62 = f - K 

63 = 11 - |n 9 

(**) = -If 

\ 102 / 

/ 19 3 44 \ 
/ 15 5 15 \ 

f 4 

I 11 I Z6 J 

\ 30 2 3 / 
/ 17 1 3 \ 

/ 10 2 2 1 

(Qa) — \ 2 2 2 

V 2 2 / 


61 = - (| + 2n g ) 
b 2 = 5 - 2n 9 
63 = 9 - 2n fl 

/__9 -9 \ 
/ 25 5 1 

{b lk )= -| 17 

\ 54 y 

/ 38 6 88 \ 
/ 15 5 15 \ 

-n g | 14 8 

U Q § 

V 15 J 3 / 
/ 26 14 18 \ 

/ 5 5 5 \ 

(Cia) =662 

V 4 4 y 
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